Low energy monopole Modes of a Trapped atomic Fermi Gas 
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We consider the low energy collective monopole modes of a 
trapped weakly interacting atomic Fermi gas in the collision- 
less regime. The spectrum is calculated for varying coupling 
strength and chemical potential. Using an effective Hamilto- 
nian, we derive analytical results that agree well with numeri- 
cal calculations in various regimes. The onset of superfluidity 
is shown to lead to effects such as the vanishing of the energy 
required to create a Cooper molecule at a critical coupling 
strength and to the emergence of pair vibration excitations. 
Our analysis suggests ways to experimentally detect the pres- 
ence of the superfluid phase in trapped atomic Fermi gases. 

Pacs Numbers: 67., 05.30.Fk, 32.80.Pj 

The trapping and cooling of fermionic atoms is attract- 
ing increasing attention in the field of ultra cold atomic 
gases. Temperatures as low as T ^ 0.2Tf with Tp denot- 
ing the Fermi temperature have recently been observed 
experimentally for such gases [1] . One of the intriguing 
goals of the impressive experimental effort is to observe 
the theoretically predicted transition to a superfluid state 
below a critical temperature Tc when the effective atom- 
atom interaction is attractive [2] . 

In this paper, the T = low energy collective monopole 
modes of a weakly interacting fermi gas in the collision- 
less limit are examined both analytically and numerically. 
We identify spectral features characteristic of the onset 
of superfluidity. Since the collective modes of trapped 
gases can be measured with high precision [3], this leads 
us to suggest appealing ways of experimentally observing 
the predicted phase-transition. 

Consider a dilute gas of fermionic atoms of mass m 
with an equal number of atoms N„ in two hyperfine 
states |<7 =t, i) trapped in a spherically symmetric har- 



monic potential of the form Uo{r) 



-/2. Assum- 



ing that the gas is in the dilute regime kp\a\ <C 1 with 
kp = y/2m^F ifJ-F is the chemical potential), the interac- 
tion between atoms in the two different hyperfine states 
can be described as gS{r) where g = Ana/rii and a is 
the s-wave scattering length. We put h ~ I throughout 
the paper. For a < the effective atom-atom interac- 
tion is attractive and the gas is unstable toward the for- 
mation of Cooper pairs below a certain critical temper- 
ature Tc [2] . The quasiparticle excitations of the super- 
fluid state can be calculated by solving the Bogoliubov-de 
Gennes (BdG) equations [4]. 

The collective excitations of the gas can be calculated 
in the coUisionless limit using the random-phase approx- 
imation (RPA). In a real space formalism, the poles of 
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yield the modes of the system. Here, {{AB)) is the 
Fourier transform of the retarded function —iQ{t — 
t'){[A{r,t), B{r' ,t')]) and (...) is the thermal average. 
Also, pa{r) = V^t (i.)^/,^(r) and x(r) = ■(Aj, (r)Vi| (r). Cor- 
relation functions with the operator x(r) are included 
in order to describe the effects of superfluidity. A nu- 
merical scheme for calculating n(cj) within the RPA was 
described in Ref. [5]. 

The nature of the superfluid state depends on two 
parameters: the number of particles trapped and the 
interaction strength [6]. When the atom-atom interac- 
tion is sufficiently weak and/or there are not too many 
particles trapped, the gas is in the so-called intrashell 
regime; it can in many ways best be regarded as a gi- 
ant nucleus. In this regime, the Cooper pairs are formed 
only between particles with angular quantum numbers 
{l,m) and {I, —m) residing in the same harmonic oscilla- 
tor shell with radial quantum number n [i.e. single par- 
ticle energy e„ = (n -I- 3/2)tJT]- For stronger interac- 
tions and/or more particles trapped such that £,bcs ^ R 
where S,bcs — kp{7TmA)~^ is the coherence length (A is 
the pairing gap) and R is the size of the system given by 
the Thomas-fermi result R ~ yji^p/m, the superfluid 
state is best described as a quasi-homogenous bulk sys- 
tem. The transition between the two qualitatively differ- 
ent regimes is roughly determined by S,bcs — Rot equiv- 
alently A ~ ojt- Using the bulk theory prediction A = 
(2/e)^/^/iF exp(— 7r/2fci?|a|) [8] and the Thomas- Fermi re- 
sult i-ip = {<(!>NcY/^ujt (neglecting small effects due to the 
pairing as A/fip <^ 1 in the dilute regime), the equation 

A = ut yields N; ^ ^ (f)"^ exp(37r/2fcj?|a|). For a given 
coupling strength with <C iV* atoms trapped, the gas 
is in the in the intrashell regime and for N^^ 3> N* atoms 
trapped it is in the bulk regime. With kp\a\ = 0.1 this 
condition yields N* = 4x 10^° and for fcF|a| = 0.3 we get 
iV* = 9 X 10^. Current experiments have N ~ 0(10^) 
particles trapped [1]. Also, the present experimental ef- 
fort uses optical traps as the conflning potential as they 
allow greater greater flexibility with respect to which 
atomic hyperfine states can be trapped [1]. The optical 
traps used so far have 1.2kHZ< WT/27r < 12kHZ corre- 
sponding to temperatures 58nK< T < 580nK [1]. Thus, 
the condition ksTc^ ojt is less restrictive than for mag- 
netic traps. Note also that fcsTc ~ A(T = 0)/2 in the 
intrashell regime can be orders of magnitude higher than 



1 



the bulk prediction [9] . The above estimates seem to in- 
dicate, that for the gas in the dilute limit it is not unlikely 
that future experiments on trapped superfluid fermionic 
atoms will be in the intrashell regime. 

For A;_F|a| > 0.3 the gas is not in the dilute limit. Re- 
cently, a number of papers have been dealing with the 
case of a Feshbach resonance mediated interaction, which 
is predicted to increase Tc dramatically [10]. This could 
make the experimental observation of the superfluid tran- 
sition easier. It is however presently unclear what hap- 
pens in this strongly correlated regime where issues such 
as the instability toward spinodal decomposition have to 
be addressed. 

The collective modes also depends qualitatively on 
which regime the gas is in. In the bulk regime, the modes 
can be calculated by looking at phase fluctuations of the 
pairing field. A gradient and frequency expansion of the 
lower right 2x2 part of Eq. (1) then yields a hydrody- 
namic spectrum for the lowest modes of the gas [7] . 

In light of the possible experimental relevance, we con- 
sider in this paper the collective modes in the intrashell 
regime. From the fact that there is only pairing between 
particles with quantum numbers {n, I, m) and (n, /, — m) 
follows a crucial simplification: Ignoring for the time be- 
ing the effect of the Hartree field, the monopole pairing 
correlations of the particles are accurately described by 
the effective Hamiltonian 



{BCS\%i.rnnnlrn^H^lvrn'^ll,V-rn'i\BCS) 
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{BCS\%i -mi'lnlm'\^n'l' -m' i'ln'l'm''\H\BCS) . The SU- 

perfiuid groimd state \BCS) and the quasi particle cre- 
ation operator with quantum numbers (n, I, m, a) 
can both can be obtained from a solution to the BdG 
equations [4]. Also, X^im = {Oljni-TnnnimiW) and 
Ynim = {Olliim^lii-miW)- t^LC intrashcU regime, 

7nim<T = UnQnlma ± l^na|j(_„_^(-l)", = {1 + 

^n/Er.)/2 and vl = {1 - ^n/En)/2 with En = (^^ + 

/\2-ji/2 j-gj -y^g have written the RPA equations appro- 
priate for subspace spanned by monopole pairing excita- 
tions of the kind + cc. 

We first discuss the solutions to Eq. (3) when 
the gas is in the normal phase (A = 0). Exci- 
tations of monopole symmetry are given by Fj, = 

Kim = (-1)™^. and Y^i^ = {-irv- in Eq. (3). Us- 
ing this we obtain after some algebra 
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in the sense that it can be shown that Eq. (2) repro- 
duces the gap equation and thus the results concerning 
the quasiparticlc properties presented in Ref. [6] in the 
intrashell regime. Here G = 32TT~'^kF\a\/15 is the effec- 
tive coupling strength and r2„y = {np + + 2)/2 
is the pair degeneracy of the highest harmonic oscilla- 
tor level occupied with radial quantum number np. The 
operator a\i^^ creates a particle with quantum num- 
bers (n, I, m, a) and energy ^„ = tn — Hf- This effective 
Hamiltonian is much simpler to work with than the full 
Hamiltonian and it allows several analytic results to be 
derived. It should be emphasized that Eq. (2) can be de- 
rived from the full Hamiltonian in the intrashell regime; 
it thus describes the actual physical system accurately 
and is not merely a schematic model as we shall see. 

In order to derive analytical results, we write the RPA 
eqn. for T = in configuration space as [11] 
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Here 

the frequency of the excitation |0) \v) = rj^|0) where 
|0) and is the RPA ground and excited state respec- 
tively. The matrices A and B are given by Anim,n'i'm' = 



and 
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Equation (4) determines the lowest collective monopole 
modes of the gas in the normal phase. These modes are 
not par- 
ticle conserving. They correspond to adding/removing 
a pair of particles; i.e. one obtains for the addi- 
tion mode ft = E„>„p;m(-l)"^"ali™T"i(- 
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E„<„pim(-l)'"^naLmT%;-mi ^1*^ energy and like- 
wise for the removal mode. The lowest particle con- 
serving monopole excitation is then given by the state 
rj;r^j|0) with energy bJr + i^a- 

Consider the case — {uf + 2)ijJt corresponding to a 
symmetric level distribution around the Fermi level with 
levels below jiF (i.e. n < hf) completely filled. In this 
case, the system is in the normal phase for T = for weak 
coupling strengths such that G < Gc- Gc is obtained 
by setting A = in the gap equation reading G~^ — 
Y^^^E-^ = 0. This yields G~^ - 25/2ln[4exp(7)nF]/3 
where 7 ~ 0.577 is Eulers constant [6]. Note that this 
equation is reproduced by setting w = in Eq. (4). We 
therefore conclude that the closed shell system becomes 
unstable toward the formation of Cooper pairs when the 
energy cost of adding/removing a pair of particles goes 
to zero. Equation (4) can be solved for G — Gc ^ 0_ 
yielding 
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with e(3) = Yfr 



,-3 



1.202. This equation gives 



the frequency of the lowest particle conserving monopole 
mode as the closed shell system becomes unstable toward 
the formation of a superfluid with increasing coupling 
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strength. For G ^ 0, we can also solve Eq. (4) obtaining 

{Wr + L0a)/0JT = 2 -AG, G ^ 0. (6) 

ftft ^ 
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We have 

for G — * 0; i.e. we simply take a pair of particles in the 
highest filled shell and put them into a monopole state 

in the lowest empty shell. This obviously costs 2ujt in 
energy for G — » 0. Equation (6) can also be obtained 
by a first order perturbative calculation of the first ex- 
cited state energy [12]. Equation (4-6) give the lowest 
monopole excitation for a closed shell system in the nor- 
mal phase. 

We now solve Eq. (3) when the gas is in the superfluid 
phase (A > 0). Using Xnim = (-1)"^„ and Ynim = 
(— 1)™!^, we obtain after some rather lengthy algebra: 
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Equation (7) determines the frequencies w of the low- 
est monopole pair vibrations of the gas in the superfluid 
phase. First, we note that it follows from the gap equa- 
tion that a; = is a solution to Eq. (7). This is the 
spurious solution coming from the fact that BCS the- 
ory breaks the U{\) gauge symmetry of the pair-field or 
equivalently particle conservation. We now examine the 
two opposite cases of interest: The completely filled and 
the open shell case. 

First consider the case oi ^ip = {np + 2)ujt with G > 
Gc such that the gas is superfluid. The solutions to Eq. 
(7) split into an even solution [X^p+n = ^np-n and 
Ynp+n = Ynp-n) Corresponding to the spurious mode 
with = and an odd solution with Xnp+n = —^np-n 



and Yn 
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solution energy as 



Some algebra yields the odd 
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where we have used A = [2--V27-i^(3)-i3]i/2(G'-i _ 
1)1/2 fyj. c > Gc in the closed shell case [6]. This 
is the pair vibration mode f p|0) where the modulus |A| 
of A oscillates. Note that 2A is smaller than the low- 
est particle conserving quasiparticle excitation energy: 

2Er,p = 2^ a, + A2 > 2A {^np = -ujt/2). The 
pair vibration mode corresponds to coherent excitations 
of quasiparticle pairs of monopole symmetry such that 
the resulting energy is smaller than 2Enp; the mode is 
strongly collective. It also conserves the number of par- 
ticles and can be excited by a particle conserving probe 
of monopole symmetry. We conclude that in the case of 
a closed shell configuration, one can detect experimen- 
tally the onset of superfiuidity with increasing coupling 
strength by observing the lowest monopole excitation 



frequency of the gas vanish as given by Eq. (5). The 
system becomes increasingly unstable toward the excita- 
tion of pairs of particles into the lowest empty band as 
they can use the degeneracy of this band to maximize 
their spatial overlap and thereby increase the (absolute 
value of) interaction energy; the energy cost of forming a 
Cooper molecule simply goes to zero at the critical cou- 
pling strength. Once the system is superfluid for G > Gc, 
the lowest monopole excitation energy again increases as 
given by Eq. (8). Now the ground state of the system 
already consist of correlated pair excitations across the 
Fermi level to create Cooper pairs with maximum spatial 
overlap and the lowest excitation is the pair vibration. 

In the case of an open shell with /j,f = {np + 3/2)ut, 
the gas remains superfluid for G — * 0+ for T = 0. In 
this limit, there is only pairing in the shell at the chem- 
ical potential {n = np) [6]. Solving the RPA equations 
for this shell, we obtain the w = spurious solution 
and the pair breaking solution with energy 2Enp = 2A 
i^np = 0). Residual interactions between between the 
excitations lower the frequency slightly below 2A. To a 
good approximation however, the lowest monopole exci- 
tation energy of the system is 

uj 2A 2G 

WT ~ WT ~ l-25/2ln(e7njr)G/3' + 
where 

we have used A/ujt = G[l - 2^/^ ln(exp(7)nF)G/3]-i 
in the open shell case [6]. Superfluidity can thus also in 
the open shell case be detected by the presence of a low 
energy a; ~ 2A mode. 

We now compare our analytical results with an exact 
numerical solution of the RPA equations obtained using 
the method described in Ref. [5]. First, consider the case 
of a closed shell system [ij,p ~ {np + 2)ujt]- In fig. 1 
(a) we plot the lowest monopole excitation energy as a 
function of the coupling strength G/Gc- For G < Gc the 
excitation is Fj^r^jO) whereas for G > Gc the excitation is 
f p|0). The X 's are obtained from a numerical calculation 
of the matrix given in Eq. (1) and the lines are Eq. (5), 
(6), and (8). We have chosen np = 50 corresponding to 
~ 2 X 10^ particles trapped in each of the two hyperfine 
states. There is very good agreement between the numer- 
ical results and the analytical predictions in their regions 
of validity. This shows that effective Hamiltonian given 
by Eq. (2) accurately describes the pairing correlations 
for the physical system in the intrashell regime. As the 
coupling strength increases, the intrashell pairing ansatz 
eventually breaks down; this occurs for G/Gc ~ 1.2 for 
this set of parameters [6] . For stronger coupling, Eq. (2) 
is invalid and one needs to work with the full Hamilto- 
nian and Eq. (1). This break down is seen numerically 
as a fragmentation and disappearance of the u> ~ 2A 
pair vibration mode for G/Gc ~ 1-2. Eventually, when 
(.BCS <C i? or equivalently G ^ Gc, the lowest collective 
excitations are not the pair vibrations corresponding to 
fluctuations in |A| but the Goldstone modes correspond- 
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ing to phase fluctuations of A [5,7]. The lowest monopole 
mode is w = 2u!t in this regime. 

Let us discuss the open shell case fip = (n + 3/2)ujt 
with n_F = 50. In fig. 1 (b), we plot the lowest monopole 
excitation energy as a function of the coupling strength 
G/u)T- The x's are obtained from a numerical calcula- 
tion of the matrix given in Eq. (1) and the line is Eq. (9). 
We see that there is good agreement between the nu- 
merical results and Eq. (9) for G ^ 0. This shows that 
in the open shell case, the lowest excitations are essen- 
tially quasi-particle excitations with energy 2A in the 
intrashell regime. The exact mode frequency is some- 
what lower than 2A as expected reflecting the residual 
interaction between the quasiparticles. With increasing 
coupling, the intrashell eventually ansatz breaks down 
for GIuINf) ~ 0.4 [6] making Eq. (2) and thus Eq. (9) 
invalid. Again, for stronger coupling when ^bcs ^ H 
the Goldstone modes are the lowest excitations. 

Several other methods to observe the superfluid tran- 
sition for atomic Fermi gases have been suggested in the 
literature. They include laser probing of the pairing en- 
ergy, looking at the expansion properties of the gas, and 
probing the density profile [13]. 

We have up till now ignored the effects of the Hartree 
field. This field introduces a dispersion on the Harmonic 
oscillator levels with respect to the angular momentum 
I lifting the shell degeneracy fin = {n + l)(n + 2)/2. 
One should therefore replace ^„ with an ^-dependent 
in Eq. (2), were S^ni includes the effect of the Hartree 
field on the single particle energies. The resulting ef- 
fective Hamiltonian describes the pairing correlations in 
the intrashell regime including the effects of the Hartree 
field. When the splitting introduced by the Hartree field 
is small compared with the pairing energy A, it can be 
ignored. This is the case for systems with < 30 (de- 
pending on the coupling strength) [6]. For very large 
systems with np ^ 200, there is only pairing within the 
levels with = and to = — — / + 1 . . . ,1. An analysis 
completely similar to the one given above yields the low- 
est excitation energy lij ~ 2A„p/ in this limit with A„^/ 
denoting the pairing energy in the level = 0. We 
plot in fig. 1 (b) as -|-'s the results for lowest monopole 
excitation energy obtained from a numerical calculation 
based on Eq. (1) including the effects of the Hartree field. 
For comparison with the results obtained neglecting the 
Hartree field, the number of particles is fixed in both 
cases. We see that for G — > 0, the Hartree field sup- 
presses the pairing energy as it introduces a splitting of 
the levels in the shell at the chemical potential such that 
the density of states for = is reduced. For stronger 
coupling however, the Hartree field increases the pairing 
energy. This is because the Hartree field increases the 
density and thus the spatial overlap between particles 
in the center of trap [6]. From fig. 1 (b), we see that 
Eq. (9) reproduces qualitatively the energy of the lowest 
mode as a function of G even when the Hartree field is 
included. This is as expected since tif = 50 is an inter- 



mediate value where the pairing is not in the limit where 
the Hartree field can be ignored, but also far away from 
the regime, where the Hartree splitting is much larger 
than the pairing energy. 

It should be noted that the conclusion that there are 
pair vibration modes with an energy ~ 2A is independent 
of the symmetry of the trapping potential. This result 
simply reflects the energy required to break individual 
Cooper pairs. However, the exact value of A as a function 
of the coupling strength depends on the symmetry of the 
potential and Eq. (5-6), (8-9) are of course only valid for 
a spherical trap. 

In summary, we examined in detail the low energy 
monopole modes of a trapped fermi g ciS clS 0j function 
of the coupling strength and the chemical potential. Us- 
ing an effective Hamiltonian describing the pair correla- 
tions in the intrashell regime we derived analytical re- 
sults that agree well with numerical calculations. As the 
gas approaches the superfluid instability with increasing 
coupling strength for the closed shell case, the lowest ex- 
citation energy corresponding to the creation of a Cooper 
molecule was shown to vanish. In the superfluid phase, 
the lowest excitation energy is ~ 2A independent of 
the position of the chemical potential. The results pre- 
sented should be relevant for the present experimental 
effort to observe the superfluid transition. 

We acknowledge valuable discussions with B. R. Mot- 
telson. 
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FIG. 1. The lowest monopole excitation energy in units of 
ljt for (a): a closed shell system as a function of G/Gc and 
(b): an open shell system as a function of G/uot- The inset 
in (a) shows the G ^ region in detail. 
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